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We exploit a novel approximation scheme to obtain a new and compact formula for the parameters 
underlying coherent-state control of the evolution of a pair of entangled two-level systems. It is 
appropriate for long times and for relatively strong external quantum control via coherent state 
irradiation. We take account of both discrete-state and continuous-variable degrees of freedom. The 
formula predicts the relative heights of entanglement revivals and their timing and duration. 
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I. INTRODUCTION 

Almost 75 years after its specific introduction by 
Schrodinger |l| a general understanding of entanglement 
within quantum theory is still being sought in the sense 
that measures of it are still being developed. Only two- 
particle entanglement is well quantified, and even in that 
case the dimension of the state space that can be taken 
into account is extremely small unless the joint state is 
pure. The behavior of few-state entangled quantum sys- 
tems takes a wide variety of forms in different contexts, 
such as cavity QED, circuit QED, spontaneous paramet- 
ric down conversion, photonic lattices and ion traps, in- 
volving atoms, ions, photons, superconducting circuits, 
quantum dots, and spins. In any of these cases, the 
controlled manipulation of entanglement with external 
agents is one approach to a form of quantum control, 
and it remains an open challenge. 

We will examine preservation rather than generation 
of entanglement, particularly the evolution of two re- 
mote qubits that are entangled and being "stored" for 
later use, ideally isolated from environmental influence, 
but subject to localized control interactions. The qubits 
have the character of a primitive quantum memory unit. 
Our approach differs further from that of most previous 
studies by focusing on the effects of (i) strong control 
fields, (ii) that are treated quantum mechanically but 
are not themselves initially entangled, and (iii) which act 
for times significantly longer than Rabi oscillation time 
scales. We will work under the condition that the mem- 
bers of the entangled pair are truly remote, i.e., do not 
interact with each other directly, or through third parties 
even indirectly. We will take the control mechanism as 
a continuous variable (CV) interaction, which is treated 
using the Jaynes-Cummings (JC) model Hamiltonian 
A sketch of a cavity QED context for this study is shown 
in Fig. [S 

There are previous theoretical discussions that also use 
the JC interaction for related considerations arising in the 
interaction of a qubit pair with CV fields. These are typ- 
ically interested in questions of entanglement generation 
in the qubits by transfer from initially weakly squeezed 
states of CV cavity fields. An early study by Kim, et al. 
3, considered the action of a two- mode squeezed field 



that is injected onto two qubits held in two separate high- 
Q cavities. From the two-qubit reduced density matrix, 
obtained by tracing over the fields, they calculated subse- 
quent atom-atom entanglement, and found that it is max- 
imized when the squeezed field state is pure. Their results 
include one of the first examples of the time-dependent 
entanglement birth and death effect (see overview [4|]). 
The question whether and to what extent the entangle- 
ment of weak squeezing could be transferred to the ini- 
tially unentangled qubits was answered positively. 

Later studies 5-7] showed that the field entanglement 
in the JC cavity could in principle be employed for trans- 
fer into networks or to subsequent pairs of atoms transit- 
ting the cavities. This work was extended [8[ to address 
an example of the "storage" problem, and transfer of en- 
tanglement from discrete qubits to an entangled coherent 
state was also examined by Zhou and Wang [9j . We com- 
ment below on relations to our results. The evolution of 
entanglement in a qubit-field system, where the qubit 
and the field start from mixed states was examined by 
Rendell and Rajagopal 1Q]. Because of the lack of an 
entanglement measure for 2 x oo systems, their calcula- 
tion of entanglement for the entire system was limited to 
a lower bound for the concurrence. 





Entanglement 



Atom A + Cavity a 




Atom B + Cavity b 



FIG. 1: (Color online) A diagram illustrating our system. 
The atoms A and B are placed in their respective cavities a 
and b after being entangled. Only A-a and B-b interactions 
take place thereafter, but non-interacting A-B entanglement 
changes in time. 



2 



We first quickly review entanglement dynamics in the 
JC context in order to establish notation for control pa- 
rameters, time scales, etc., and to allow direct interpre- 
tation of a new approximate formula - see relation (|60[) 
below - for entanglement control parameters. We restrict 
attention to the dynamics of entanglement between two 
sites holding qubits A and B. We suppose that there is 
quantum memory present in the qubits that is of interest, 
i.e., that the two atoms have been entangled in some way 
before being inserted into their respective cavities. Field 
modes a and b are available at the A-B sites that can be 
used for interacting with the atoms (i.e., for "controlling" 
them and their entanglement externally). This scenario 
is suggested in Fig. [1] 

II. HAMILTONIAN EQUATIONS 

The Jaynes-Cummings |1| Hamiltonian (with h = 1) is 
given by: 

tt ^0 A f # + A A \ \ 

Htot = -^-cf z + (g A a< cr_ + g A a + a) + Lja'a 

+Y°? + {9* B b^ B +gB^b)+^b, (1) 

where ujq is the transition frequency between the two lev- 
els of the atoms, g is the constant of coupling between the 
atoms and the fields and uj is the angular frequency of the 
single-mode field. The usual Pauli matrices <rf' B , <r+' B 
and a^' B describe the atoms, while a\ a and b\ b are 
the raising and lowering operators for the fields. The dif- 
ference between g and g* in both cases is only a phase, 
but there may be a substantial difference between g A and 
gs if the cavities are different. 

The Hamiltonian does not support any interaction be- 
tween atom A and atom B or between mode a and mode 
b. Physical realization of our scenario does not appear 
out of the question, as the Jaynes-Cummings model has 
been realized in the laboratory in several well-known 
ways 

To illustrate the approach to be followed when the 
fields are taken as highly excited (nearly classical) co- 
herent states, we start with a simpler example, in which 
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FIG. 2: (Color online) Sketch indicating non-interacting 
qubits in a quantum storage network. Quantum memory in 
the form of two-qubit entanglement between any two rela- 
tively remote sites, indicated here by a dashed line, can be 
altered by fields local to the sites. 



the field modes are initially in their vacuum states: 
|0 o ) ® |0h), and the two atoms are in a superposition of 
Bell states. An additional simplification, to focus on the 
most elementary aspects of the entanglement evolution 
first, will be to take the cavities identical, and assume 
9 A = gs = <?• For the vacuum case the coupling terms 
in our Hamiltonian indicate that for creation of a pho- 
ton in a cavity the atom in that cavity must decay to 
a lower state. Since the cavities each contain only one 
atom, there can be no more than one photon in each at 
any time. This means that the cavities are also two-state 
quantum systems, i.e., qubits. 

The eigenstates of the JC Hamiltonian are well known. 
We will write for either Aa or Bb 

Hjc^t) = A±h£>. (2) 

We will denote the excited and ground atomic states by 
|e) and \g), and denote the cavity modes' photon states 
by the photon number n. Then the JC eigenvalues are 
given by: 

A±=nw + i(A±VA 2 + G»), (3) 

where A = ujq — lu is the detuning and 

G n = 2g s/K, (4) 

is the n-photon Rabi frequency. 

The JC eigenstates have the following form as super- 
positions of the bare atom and cavity states: 

l^o) = Is; 0) (5) 
= c n \e;n- 1) +s n \g;n) (n > 0) (6) 
IVO = -s n \e;n - 1) + \c n \g;n) (n > 0). (7) 

In these equations we have introduced some convenient 
abbreviations: 

c n = cos(6>„/2) and s n = sin(0„/2), (8) 

where the rotation angle 6 n can be identified with the 
Bloch sphere polar angle and is defined in the usual way: 

cos# n = — z and sin6>„ = — = =. (9) 

In this preliminary example we will need the true 
ground state |^o) = |<?j 0) and the two dressed states for 
n=l. These three states are closed under the JC Hamil- 
tonian for each site. In other words, we use only n = 1 
in the equations above, so the subscript n can mostly be 
ignored and we will frequently drop it (A n — > A, c n — > c, 
etc.). 

III. MEASURE OF PAIRWISE 
ENTANGLEMENT 

The JC dressed states are atom-photon locally- 
entangled states themselves, and this entanglement has 
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interesting consequences if the cavities are fairly highly 
excited, as Gea-Banacloche [l5[ and Phoenix and Knight 
16] originally discussed. Here we are interested in non- 
local atom-atom entanglement between the two spatially 
separated sites on a qubit lattice, as suggested in Fig. [2] 

In the general context of entanglement we note that 
there is no accepted and practically workable criterion for 
determining separability of arbitrary four-particle states. 
Our purposes will be satisfied by working with the two- 
qubit atomic states obtained from the time-evolving four- 
qubit state. All familiar measures agree about separabil- 
ity in the two-qubit domain of entanglement. That is, 
entropy of formation, Schmidt number, tangle, negativ- 
ity, and concurrence are numerically somewhat different, 
but in the two-qubit domains of their applicability they 
are in full agreement when they signal entanglement or 
lack of entanglement. 

We adopt Wootters' concurrence [ItJ as our measure 
in this discussion, mainly for its convenient normaliza- 
tion: 1 > C > 0, where C = indicates separability 
(zero entanglement) and C = 1 means maximal pure 
state entanglement, as in a Bell state; and simplicity of 
calculation: 



C(p) = max{0, v Ai 



A 3 - VA 4 }, (10) 



where the quantities A^ are the eigenvalues in decreasing 
order of the auxilliary matrix 



C = P((7y <g) (Ty)p*(a y <g) (Ty), 



(11) 



where p* denotes the complex conjugation of p in the 
standard basis and a y is the Pauli matrix expressed in 
the standard basis. 

One finds that in reduction to two-qubit form, by trac- 
ing over the two field-mode qubits, the resulting two- 
qubit mixed state always has the X form [l8( , where only 
the diagonal and anti-diagonal elements are not zero: 



pn P14 

P22 P23 

P* 23 P33 

PlA p 44 



(12) 



where pn + P22 + P33 + Paa = 1- The concurrence of this 
mixed state is easily found to be 

C = 2max{0, |p 23 | - ^pnpa, \pu\ ~ v 'P22/O33} 
= 2max{0,Q}, (13) 

so it is clear that Q, which is the larger of \p23 \ — \fp\iP~4A 
and I pi«4 1 — */ 'P22P33, will be an important quantity. 
For example, Q{t) obeys certain conservation relations 
(whereas C does not) in some special cases because it can 
be negative. Furthermore, one can see [l9[ that while 
Q < implies a separable state (C = 0), the slightly 
stronger condition Q < implies both C = and also 
that the state is mixed rather than pure, information not 
available from C. 



IV. PARTIALLY ENTANGLED BELL STATES 

In the two-site situation under consideration there are, 
in principle, six clearly distinguishable entanglements. 
They have been overlooked in the earlier papers men- 
tioned @,[1-G3, but they all can carry information about 
the bipartite entanglements that may arise. With an 
obvious notation we can denote these concurrences as 
C AB , C ab , C Aa , C Bb , C A \ C Ba . Except for C Aa and 
C Bb , they measure remote (non-local) entanglements, 
and their dynamics have been reported to be restricted by 
previously unknown invariant combinations of entangle- 
ment parameters (20l - [22l | . For background we summarize 
results when the initial states are superpositions of the 
Bell states: 1$^) ~ \e A ,e B ) ± \9At9b), which we write: 



|$Ab) = cosa\e A ,e B ) + sma\g A ,g B ) 



(14) 



It is easy to see that a = ±7r/4 reproduces the two pure 
Bell states, and other values of a represent an adjustable 
phase between the entangled atomic states. The optional 
value of a is significant because it can be understood as 
a measure of the relative difference in preparation of the 
two atoms in the two identicial cavities. This can be im- 
portant from the point of view of experimental relevance, 
as identical preparation of the atoms would be difficult. 
The initial state for the atoms plus cavities is therefore: 

|*(0)> = |0 a ,0 b ) (15) 

= (cosa|eA,e B ) + shia\g A ,g B )) <g) |0 a ,0 6 ). 

In terms of the dressed eigenstates given above ([5]), we 
can rewrite: 

|e A ,0„> = c|Vf) - «|Vf) 

\g A , l a ) = + <#i~) and 

\g A ,0 a ) = |Vo). (16) 

Thus the initial atom-atom entangled state has the form 

cosa|eA,0 a ) ® \e B ,0b) +sma\g A ,O a ) ® \g B ,0b) 
= cosa(c\^) A - s\i/)i) A ) (<#i")s - s\ipi) B ) 
+ sma\iP°) A ® \ij)°) B . (17) 

Evolution in time is easily followed through the evolu- 
tion of the dressed states: 



\i> ± {t))=e 



|^(0)). 



(18) 



Note that since the combination of coefficients in |$(0)) 
uniquely associates c with \ip + ), and s with the time 
evolution can be transferred to the c and s symbols. We 
will henceforth consider them carrying the time-evolution 
exponents. We will use the notation cq = c(0) and sq = 
s(0) to refer to their values at t = (no relation to the 
n = subscripts in Eq. ©), so that, 



c = c(t) = c Q e 



s(t) 



s e 



i\+t 



-i\~t 



(19) 
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where A + and A~ are obtained by inserting n = 1 into 
Eq. ([3]). Then we can write (temporarily again indicating 
explicit time dependences for the c's and s's): 

|*(t)> = cosa(c(t)\i;+) A -s(t)\^) A ) 

<8 (c(t)|V»f>B -s(<)IV'Ds) 

+ sina|i/; ) A ® |^ > B , (20) 

where |'0 ± ) will continue to refer to the JC states at t = 0. 

Now we revert to the bare basis states in preparation 
for the tracing needed to calculate Q AB and C AB : 



|$(t)) = cosa\c(t)(c \e A ,0 a ) + s \g A , 1 )) 

- s(i)(-s o |eA,0 a + coIsa, la))) ® fc(f)(co|eB,0 6 

+ solas, Ifc)) - s(t)(-s |eB,0 6 + c Q \g B , h)) 
+ sina\g A ,0 a ) ® |.g B ,0 o ) 
= cosa((c(t)c + s(t)s )|e J i,0 o ) 

+ (c(f)so - s(t)co)|5A, 1 Q 
® Uc(t)c + s(t)s )|e B ,0 o 

+ (c(t)so - s(t)co)|flfl, lfe 

+ sina| 3A ,0 a ) ® |.g B ,0 o ). (21) 

To get the two-qubit mixed state needed for calculation 
of Q AB the projections that are needed are: 

(0 a ,0&|$(t)) = cos a(cc + ss ) 2 \e a , e B ) + sin a\g A ,g B ) 

(l a ,0 6 |$(i)) = cos a(cs - sc )(cc + ss )\g A ,e B ) 

(Q a , 1&|$(*)) = cosa(cc + ss )(cs - sc )|eA,5s) 

(l«,l&|*(t)) = 0, (22) 

where now c and s have replaced c(t) and s(t) to conserve 
space. These results show that the AB mixed state has 
the X form mentioned in the previous Section: 



pu pu 

w; p 2 2 

p 33 

P*4 (044 



for which the concurrence has the stated form 
C AB = 2 max{0, |pu| - Vp^i}. 
and we easily find the following 



(23) 



(24) 



| Pi4 1 = | sin a cos a | (c + s + 2cqSq cos 5t), 

b = c = cos 2 a\cco + sso\ 2 \csq — sco\ 2 

= cos 2 a(cg + Sq + 2cqSq cos 5t) 

x CqSq(2 — 2cosi5i), 



(25) 



For simplicity we will evaluate this in the resonance 
case, 6 n = w/2, where c = s = l/s/2 and A = 0. Then 
we find 

1 



\pu\-y/P22p33 = t cos a[2+2cos(5t)][| tana|-sin (gt)], 

(26) 

from which the expression for concurrence is found to be: 



C 



AB 



2max{0,Q A±f }, 



(27) 



where Q AB = cos 2 a cos 2 (gt)[\ tan a\ — sin 2 (gt)}. 

Fig. [3]shows that for aj^ 7r/4 the C AB curves have the 
"sudden death" feature That is, the entanglement 
non-smoothly becomes zero and stays zero for a finite 
interval of time. The death and rebirth dynamics of non- 
interacting entanglement shown is however still governed 
by the same Rabi time scale that one associates with 
vacuum Rabi oscillations [23, urn . 




where 5 = \/ A 2 + 4g 2 . 



FIG. 3: Time dependences of AB and ab concurrences for 
three values of the superposition angle a. Note that almost 
all curves show a time interval over which C=0, i.e., during 
which the underlying state must be separable. The time scale 
is the vacuum Rabi period. From ref. [2p| ]. 



V. OPEN-SYSTEM TWO-QUBIT THEORY 

Realistic quantum control almost necessarily implies 
engagement of continuous variables and the interaction 
of qubits or other systems having a finite number of 
states with one or more "large" systems with continu- 
ous quantum states. We can use the JC formalism de- 
scribed above to enter this domain by introducing coher- 
ent state fields at the AB network sites. In Fig. 2] we 
indicate a post-selection method for experimentally ob- 
taining a Bell State to work with, in principle explaining 
the entanglement stage left open in Fig. [1] That is, the 
cavity shown in Fig. 0] is used to prepare atoms to be 
inserted into the two cavities in Fig. [TJ It is initially 
prepared in a single photon state, e.g., by micromaser 
methods. After the atoms pass through it, each entering 
one of the final cavities shown in Fig. [T] to begin their 
JC interactions, this cavity is monitored for several life- 
times, to ensure that its photon was transferred to one 
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FIG. 4: (Color Online) Top and bottom sketches show a pair 
of atoms entering and then leaving a prepared resonant cavity, 
indicating in principle a preparation of the Bell State \P ab in 
advance of insertion into the two cavities that are shown in 
Fig. [T] This is a post-selection approach, as described in the 
text. 



of the atoms A and B (without knowledge which one), 
and hasn't remained in the cavity. That run of the ex- 
periment is cancelled if /when a cavity photon leaks out, 
and the experiment must be restarted with another pair 
of ground-state atoms. 

We retain almost all of the simplifying approximations 
made earlier, and add one more by taking the two fields, 
now modeled as coherent state fields, to have the same 
average photon number n, which could be arranged ex- 
perimentally by feeding both cavities from the same co- 
herent state laser field via a 50-50 beam splitter. How- 
ever, one important simplification that we relied on pre- 
viously must now be discarded. The coherent-state fields 
have many occupied photon- number states, so the cavi- 
ties will not be two-state qubits. We assume initial atom 
entanglement in the form of the pure Bell State prepared, 
e.g., as in Fig. [4J 

\9AB{0)) = (\eg) + \ge))/V2, (28) 

and write the field state as the coherent state product 

I*o6(0)> = \a) (8 \a). (29) 

As a result our initial state for the whole system is, 

|#tot(0)> = |*as(0)> ® \a) ® \a). (30) 

The coherent states are given by the familiar Fock state 
expansion 

oo oo _ | Q | 2 /2 n 

\a)=J2 A n \n) = e n ° In)- (31) 

— n — n VU 

n— n— v 

Now the photonic density matrix is infinite- 
dimensional, and the joint AB dynamics extremely 
complicated. Fully numerical analysis is possible, but 
the insights from analytic results are highly desirable. 
We have found a key step permitting analytic progress. 
This is an apparently drastic simplification of the con- 
tinuous state spaces of the two field modes. We assume 



that it is satisfactory to replace \a) by \n). This Ansatz 
is at least weakly supported by the knowledge that 
photon number in a coherent state is Poisson-distributed 
and relatively tightly centered around n when n 3> I. 
Thus we represent the initial field state by the single 
product Fock state \n) ® \n). Note that during the JC 
interaction the photon number in either field mode a or 
b can then be n or h + 1 or n — 1 . 

Under the Ansatz mentioned, by tracing the field mode 
states, we find that the reduced density matrix for the 
qubits becomes: 

/ pxi x x x \ / pn \ 

I X P22 P23 X I _^ I p 2 2 P23 | 
I X p 2 3* P33 X I I P23* P33 I ' 

\ I X X P44 ) \ /?44 / 

where we have used the standard two-qubit basis 
[ee, eg, ge, gg}. The elements indicated by x are zero 
because of the equal-n simplification. Thus, under the 
assumptions mentioned, p is again of X-type [HI). 
For this A-type p, concurrence turns into: 

C{p) = 2 max[ 0, |p23 1 - s/pnPu }■ (33) 

The coherent fields induce time-dependent change of the 
elements pn and P44, and their growth and any decline 
of P23 will cause entanglement to decrease. 

In order to calculate the time evolution of this state 
under the JC Hamiltonian we need to calculate the time 
evolution of the the states of the individual sites, i.e., 
|e) ® \n) and \g) ® \n) for all n. The time evolution of 
these states for site A (and similarly for site B) is given 
by, 

e iHlt \e;n) = cos(gWn+ I)|e;n) 

- isin(gt\/n+ l)\g;n+ 1) (34) 

e tH ^\g;n) = cos(gt^/n)\g;n) 

— i sm(gt\/n)\e; n — I). (35) 

Using these results, the time evolution of the initial 
state of the system is found to be, 

|*tot(t)> = e^ 4 |* tot (0)> (36) 
= —7= ^ ] ^ ] A n A m ^K mn ^J, (37) 

V* n= O m= Q 

where 

K mn = —iC n+ iS m \e,e;n,m—l) + C n +iS m \e,g;n,m) 

- S n+ iS n \g,e;n+ l,m) - iS n+ \C n \g, g; n + l,m) 

- iS n C m+ i\e,e;n-l,m+l) 

- S n S m+ i\e,g;n- l,m + 1) 
+ C n C m+ i\g,e;n,m+ 1) 

- iC n S m+ i\g,g;n,m + l), (38) 
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where C n — cos(gt^/n) and S n — sm(gt^/n). The density 
matrix for the system is then given by 



Plot = \^tot{t)){\^tot{t)\. 



(39) 



and the reduced density matrix for the atoms, pab, is 
given by p AB = Tr^ m )p tot . 



VI. FORMULAS FOR CONTROL 
PARAMETERS 

It is obvious that the quantities denoted P23, Pn, Paa 
are the control parameters for atom entanglement. Hav- 
ing used the Fock state shortcut to obtain (|32]). we now 
avoid using it further and introduce the approximation 
method that leads to analytic formulas for P23, Pn, Pa 
appropriate to pab for the coherent state. 

One can show that the p 2 3 term is given by the doubly 
infinite summation, 

n,m 

— A n A n —iA rn A m j r \S n C n j r \C m S m j r \ 
+ A n A n -2A m A m+ 2S n S n -iS m +iS m +2 

— A n A n -iA m A rn+ iS n C n -iS rn -i-iC m +2 \, (40) 

Similarly the series summations for pn and pa are: 



a — 2 { A^lCn+l^n 



A n A n + 1 A m A m _ i S n 4. i C n 4. i S m C m 

a2 a2 q2(~<2 

A n A n ^ iA m A m +iS n C n S m +iC rn -i-i \ 



and 



A2 a2 q2 



/I2 A2 c-2 fi2 



+ A n A n +iA m A m —iS n +iC n -i-iS m C m 

i /i 2 yi2 ^2 ct2 

+ A„A„_i^ m A m+1 5'„C„5 m+ iC m+ i|. (42) 

The infinite extent of these summations of course reflects 
the fact that we are dealing with a quantum open system, 
by having coupled the qubits to an infinite state space. 

The sums above cannot be completed, but excellent 
analytic approximations can be found for coherent states 
that are even only moderately strong, i.e., a > 10. We 
will use the familiar Stirling formula for n!, 



i! = V2Tmn n e~ 



(43) 



and Euler's formula to approximate the terms in the sum- 
mations above by integrals. We begin by approximating 



the terms like A n A n+ iA m A m _i with A 2 n A^ which in- 
troduces an error of order 1/n near the Poisson peaks 
n w m ~ n. Then we get for z 



P23 



n 

— 2^ AnA^S n C n +iC m S m +i J 

n 

In the same way, Eq.(|4T|) and Eq.(|42|) become, 



(44) 



(45) 



]T A n C 2 n ) ( A 'X) + ( J2 AlS n C n ) .(46) 



P44 



Note that Eq. (g5]) and Eq. flU]) imply that p u = p 44 
whenever our large n approximation is valid. Now we 
rewrite C n C n +\ as, 

C n C n+ i = i| cos[gt(^/n + \Jn + 1)] 

+ cos[5*(Vn + 1 - Vn)]}, (47) 

and use the peaked nature of A n to focus on those terms 
near to n to introduce the approximation 



y/n+1 = + 



(41) which compresses C n Cn+i to 



CnC„+i — — 

Similarly, 

S n C n +l — 
S n +lC n = 



cos(2gty/n) + cos 



2v^ 



(48) 



(49) 



l r 

2 
1 
2 

1 r 



gt 

2Vn 



sm(2gt^/n) — sin 
Cm(2gt^/n) + sin ( 



cos(2gty/n) 



gt 



2y/n 

gt 

2^i 



(50) 



With these results we can simplify /?23 further: 

n v n v 

+ (jTA 2 nC os{2gt^i)) 2 

n 

-(^^sin(2 ff ^)) 2 ]. (51) 



7 



Now, using the identities 
C 



1 + cos(2gty/n) 



2 1 - cos(2gt^/n) 
2 

we can rewrite p\\ and ^44 as, 

P\\ = P44 = jjl- (^2Alcos(2gt*/n) 

n 

+ (^A2 s in(2 5 ^)) 2 . 

n 

Then Eqs. ^ and ([53]) lead to: 

P23 ~ VPHP44 = ^[(Xl A « cos (^7^)) 

n * 

5* N 2 



(52) 



(53) 



yA 2 sin(-^) 
^ " 2^ 

A 2 cos(2 5 i^) 

n 

2(^A 2 sin(2 5 V^) 



- 1 



(54) 



We can calculate the sums involved here by rewriting 
them as integrals, treating the integer n as continuous, 
again relying on the large value of n. The first two inte- 
grals we need to calculate are, 

h= A 2 . cos(-^-=) dn, and (55) 
Jo V™ 

gt 



h= I A 2 sin(-^)dn. 



(56) 



We will combine these integrals, I\ + il 2 = I12 in order 
to work with the exponential of igt/2^/n. This, together 
with Stirling's approximation on A\ , and the abbrevia- 
tion t = gt, leads to 



v 2 a 2n e" 



(57) 



The saddle point method is appropriate to calculate 
this integral, and some details are reserved for the Ap- 
pendix. The expression for I\ + il 2 is found to be, 

^exp(-^) e W- (58 ) 

Then helpful cancellations can be identified, and we 
obtain an approximate expression for |p 2 3| — \JP11Pa, 
the entanglement determiner: 

1 J e - 9 2 t 2 /Hm 2 _ 1 



\P23\ ~ VP11P44 = 7 le 



n 

i[^^ 2 sin(2gt^)] 2 . (59) 
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FIG. 5: The analytical and numerical results for concurrence 
of two qubits exposed to two quantum-coherent driving fields 
with a = 10. The evident revivals are predicted reasonably 
well by the approximate formula (|60|) . Analytical results are 
for the X-type p while the numerical ones are for the original 
P- 



The summations in (|59j) involving cos(2gt^/n) and 
sin(2giy / n) can be combined into a single exponent con- 
taining the argument 2igt^/n, which is similar to that in 
([57| , but the resulting saddle point analysis is more com- 
plicated because now y/n is in the numerator. Details are 
relegated to the Appendix, which leads to the following 
working formula: 



\P23\ - VP11P44 = ^{ CX P ( 

3(4ar)| 



16a 4 



-1 + e T /2 cos! 



y — 



exp 



2(r - 271-fca) 1 



fc=l,2, 

x cos[4a(T — 27rfca)] 



1 + 7T 2 fc 2 



(60) 



In writing this last expression we have used the fact that 
around r = 2irka only the term with the corresponding k 
gives a significant contribution to the sums. The contri- 
bution to r = 27rfca from any other k' is proportional to 

exp ( — 47r 2 a 2 (fc — k') 2 /[l + 7r 2 (fc') 2 ]^ , so it decays expo- 
nentially with the distance from k. Thus the final term 
is the main result, which can be read separately for each 
value of the step index k. 



VII. OVERVIEW AND IMPLICATIONS 

The X-state simplification and Fock-state Ansatz we 
introduced work together in a surprisingly accurate way. 
The evaluation of concurrence for n 3> 1 here is generi- 
cally the same as that presented for zero detuning qubit 
inversion in the original discussion of quantum revivals 
26] , and revivals arise in Fig. [3] because of our focus 
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on strong control and long-time dynamics. The figure 
shows both the predictions of the analytic expressions 
given above, and also the results of a numeric check of 
the Fock state Ansatz that initiated the analytic calcu- 
lations. 

While entanglement recovery is never complete, the 
question how complete is important, and it was posed 
by Lee, et al. Q. We can display the answer in another 
analytic formula, giving relative revival envelope heights: 

f 1 - exp(-T 2 /16a 4 ) 
~k 2 ■ (61) 

where k is the revival number. 

We also note that while revivals are quite robust, the 
exact-numeric to approximate-analytic agreements are 
not perfect, and the differences between them are illu- 
minating. We refer to the absence of Rabi-type behavior 
during the revivals in the numeric results, which is not 
predicted analytically. That is, the ESD events (returns 
of concurrence to zero) that occur within the analytic re- 
vival envelope are not present in the numeric envelope. 
This is shown in the expanded views of the first revival 
in Fig. |B1 

The absence of Rabi-type behavior in the numeric 
curves is seen clearly by inspecting a revival envelope in 
detail, as is done in Fig. [6] The analytical results show 
rapid Rabi-type oscillations with period r = ir/(2g^/E). 
The analytic formula retains the entanglement death and 
rebirth episodes that occur on the rapid Rabi-period 
scale, as were shown in Fig. [3] and have been discussed in 
the literature repeatedly for few-photon excitation. By 
contrast, the numeric results show a smoothed version 
without rapidly recurring ESD events. Even the zero re- 
vival, i.e., the period that is referred to as the Cummings 
collapse in the inversion literature, shows no ESD events 
within it. 

Revivals have been demonstrated experimentally for 
values of n > 3 (see, e.g., [Hj]). Our results extend the 
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FIG. 7: Similar to the bottom curve in Fig. except that 
here a = 5 and 6, rather than 10. 



revival results reported briefly by Paternostro, et al. 
which were restricted to shorter times and values of n too 
small to exhibit the wide revival separations in Fig. [5j 

Finally, the quasi-periodic modulations evident in the 
numeric details in Fig. [6] deserve comment. For contrast, 
we show in Fig. [7] below the numeric revival details for 
two other values of coherent state photon number: h = 
25 and 36. In those graphs modulations also appear, but 
with different main periods, their frequency increasing 
linearly as a increases, viz., n main modulation periods 
for a = n over a unit interval in gt/n. These modulations 
are an artifact of the assumption that no difference exists 
between the n values at the two sites being managed. 
This is an instance where a treatment with more options 
than presented here (e.g., one with different ns) should 
be expected to promote a desirable feature, the decrease 
in mid-envelope modulations. 



Analytical 
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FIG. 6: Details of the first revival shown in the preceding 
figure. The deep modulations in the analytic envelope are not 
present in the numerical envelope, which retains entanglement 
robustly through the revival event. 
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IX. APPENDIX 

Here, we sketch the saddle point analysis of the re- 
maining sums in Eqn. (I59|) . They can be combined, and 
converted to an integral: 

/>oo 

hi= / exp[a 2 f(n)]dn, (62) 
Jo 
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where 

f(n) — — (2n In a — — ln(27m) 

- nln(n) +n + 2iry/n\ - 1. (63) 

In the saddle point method, an integral of type 
jt> e Mf{x) dx ^ where M 

is a large number and f(x) is a 
twice-diffcrentiablc function, is approximated by, 



Thus, we have, 



e M ^dx 



2tt 



M\f"(x )\ 



,Mf(x ) 



(64) 



where xq is the global maximum of f(x). In our appli- 
cation Eq. (|62p can be written in this form if we choose 
M = a 2 and f(n) as given above. For the saddle point 
maximum of f(n) we need to find the point no where 
/'M =0: 

/>o) = 4f21na--^ ln(n Q ) + iTn- 1/2 ) = 0. (65) 

or V zno ) 

Assuming that \n§\ is large, we can replace ln(no) = 
ln(a 2 ) + irriQ 1 ^ 2 , and then letting n = pe l9 gives 

lnp + i6 = ln(a 2 ) + iTp- 1/2 e~ i6/2 , (66) 

where equating the real and imaginary parts of both sides 
gives, 

ln(p) = \n(a 2 )+Tp- 1 / 2 sm(9/2) 

9 = t P - 1 / 2 cos(9/2). (67) 

For r = 27rfca, where k is an integer, the equations 
above become, 



p = a 2 and 9 = (-l) fc 27rfc. 



(68) 



Now, let t — tq + At and 9 = 9q + A9, where To = 27rfca 
and 9o = 2nk. Assuming that both At and A9 are small, 
we can write 

An 

sin(0/2) = (~l) fe — and cos(<9/2) = (-l) fc . (69) 

Then Eq. (JBTJ) turns into, 

p a 2 [l + (-l) fe 7rfcA(9] 

nkA9 , _ t Ar 

Q 



= va 



+ (-l) k — . 



(70) 



In order to arrive at the second line we have used, 

,7rfcA0\ 



-1/2 



o- x (l-(-l) fc - 



2 r 



(71) 



in the second line of Eq. (|67p and we have ignored the 
terms with AtA9. 

Now, the last two terms in the second line of Eq. ([70)1 
are just A9, so, 



E*^ + (-l)^ = A*. 
2 a 



(72) 



A<9 



(-1)* 



a(l + 7r 2 fc 2 ) 
Using this in Eq. (I7D|) . we obtain 

1 



At, 



(73) 



p = a" [1 + irk 



-1) 



a(l + 7r 2 fc 2 ) 
1 



At] 
At. 



a(l + 7r 2 fc 2 ) 

Now, by inserting n = pe l8 in Eq. (|88[) . we find 



(74) 



« 2 /K) = VMp/« 2 )-Mp 1/2 )- 

pe ie (l-i0)+2vrp l ' 2 e ie ' 3 
- a 2 - ln(-\/27r). 



(75) 



The real part of this equation is given by, 

i?e{a 2 /(7i )} = -pcos9\n(p/a 2 ) - ln(p 1/2 ) + p9sm9 
+pcos9 -2rp 1/2 sin(6>/2) 
-a 2 - ln(V27r). (76) 

Using 9 = 9q + A9 and retaining only the terms upto the 
second order in A9, 



Re{a 2 f(n Q )} = -p 



1 - 



(A9f 



ln(4 

or 



Hp 1 ' 2 ) + pA9(9 + A9)+p 



-{-l) k Tp 1 ' 2 A9 -a 2 - bi(v^7r). 



{A9f 



(77) 



Writing A9 in terms of At and ignoring the terms after 
second order, 

He{QVW}a __i_(l±|^!) (AT)2 

- InGo 1 / 2 ) - Hp 1 ' 2 ) - \n{a^ 2 ) - Ih(v^tt). (78) 

The imaginary part of a 2 f(no) is, 

Im{a 2 f(n )} = - sm0p\n(p/a 2 ) - 9/2 
-p9cos9 + p sin 9 + 2rp 1/2 cos(6>/2) . (79) 

Again, using 9 — 9q + A9 and t = tq + At, writing AO in 
terms of At and ignoring the terms after second order, 



Im{a 2 f(n )} = {~l) k {2nka 2 + 



2a 



1 



2a(l + 7r 2 fc 2 ) 



irk 



1 + ir 2 k 2 



1 



2(l + 7r 2 fc 2 ) a 2 



(Ar) 2 }. 



At 

(80) 



For k — the (At) 2 part of this equation vanishes. For 
k = 1, 2, ... this part can be ignored as well. Thus we are 
left with, 

Im{a 2 f(n )} Si (-l) k (2irka 2 + 2aAr). (81) 



10 



We also need to calculate ^/27r/a 2 |/"(n )|. By using 
Eq. §B, we find 



n«o) = i( 



1/1 1 irn. 



-3/2 



a 2 V2np n 2 
which has two forms: for r = 0, 



(82) 



2n 



V « 2 I/"MI 

and for r = 27rfca (k=l,2,...): 



^ V2T 



2tt 



27TQ! 



Ya 2 |/"(n )| V ^fc 
Then the integral in Eq. (|62[) is given by, 



2vr 



'34 



a 2 |/"(n )| 



" 2 /(«o) 



(84) 



(85) 



In order to find the value of this integral we should add 
the contributions from all k = 0, 1, 2, .... As a result, the 
integral is, 



7 34 = e- T -/ 2 e 2laT 

+ \/ ~r ex P 

k=l,2,... 

x cos[2a(r — 27rfca) 



(t - 27rfca) 2 



l+7T 2 fc 2 



(86) 



Now, inserting the values of all four integrals into Eq. 

si, 



\P23\-Vpnm = ^{ ex p( 

2 v-^ 1 r 



16a 4 
exp 



e r cos(4aT) 



fe=l,2,.. 

x cos [4a (t - 2 



(- 

7rfca)] I 



- 1 

2(r - 2?rfca) 2 



1 + 7T 2 fc 2 



(87) 



Writing this last equation we have used the fact that 
around r = 2nka only the term with the corresponding 
/c gives a significant contribution to the squares of the 
sums in Eq. (1541) . The contribution to r = 2nka from 
any other k' is proportional to exp(— 47r 2 a 2 (/c — k') 2 /[l + 
ir 2 (k') 2 ]), so it decays exponentially with the distance 
from the peaks identified with integer k. Thus, while 
taking the squares in Eq. (|54[) we can ignore the cross- 
terms. 

Now, we turn to the integral given in Eq. ([57| . We 
will continue to use the saddle point method. This time 
we need to find the maximum of the function: 



f(n) = ^(2nln(a)-iln(27m) 



For this, we need to find the point no where /'(no) = 0: 

/'(n„) = ^ ( ln(a 2 )-^-ln(n )-^n " 3/2 ) =0, (89) 

Again assuming that | no \ » 1 , 

,2 



(83) Letting 



then, 



ln(n ) = ln(a 2 ) - 3/2 . 



n = pe , 



ln(p) + i9 = ln(a 2 ) - ^p^ 2 [cos(y ) - 1 sin(y )], (92) 

and by matching the real and imaginary parts of both 
sides, we find two coupled transcendental equations: 

Hp) = i n(a 2 ) _ L p -m 



(90) 
(91) 



,36> 



(93) 



We are going to retain only the terms up to second order 



m rp 



-3/2 



T , 2 , W 2 , 

--^3 and P = « 2 (l + — ) 



(94) 



With this restriction, and after inserting n = pe into 
Eq. 0, we find 



a 2 f(n ) 



-ln(a 2 ) - ln(a) - m(V27r) 
2 ^ 



(95) 



Writing p in terms of and retaining only the terms up 



to 



a 2 /(n ) 



-ln(a 2 ) - ln(a) - ln(V2^) 



W 2 a 2 9 2 


t8 


4 + 2 


4a 


f t e 


7t9 2 


+ \2a-2- 


16a 



(96) 



so that by inserting — r/4a 3 for 9, we obtain 



a /(n ) = -Zn(a J ) - ln(a) - ln(V27r) 



32a 4 



3r 2 
16a 6 



. / T T 

A2a ~ 8a 5 " 



7r 3 



256a' 



(97) 



In order to find the final form of Eq. (|64p we need to 
calculate |/"(no)| as well. 



-n ln(n) 



1. 



,/// \ 1 ( 1 1 3ir _ 

/ («o) = — ( r-f + -3-™ 

a^ \2nQ no 8 



5/2> 



(98) 



11 



Retaining the terms up to a 2 in the parentheses, we 
have |/"(n )| ^ ^. Then we use Eq. §7$ and Eq. §Z§ 
and insert a 2 for M, to get 



For t ~ a 2 we can retain the first terms in the parenthe- 
ses and ignore the rest since a 2 >> 1. Thus, we finally 
find 



2vr 



« 2 |/"M| 



e a S exp 



32a 4 



exp 



V2a 



IT 
8^ 



3r 2 
16a 6 

7ir 3 
~ 256a 7 



9) 



Zi2 = exp 



(~3^?)' 



(100) 
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